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Abstract 

We showed that if the non-thermal emission from the Galactic center in the range 
14 — 40 keV is due to inverse bremsstrahlung emission of subrelativistic protons, 
their interactions with hot and cold fractions of the interstellar medium are equally 
important. Our estimation show that about 30% of the total non-thermal flux from 
the GC in the range 14 — 40 keV is generated in regions of cold gas while the rest 
is produced by proton interaction with hot plasma. From the spatial distribution of 
6.7 keV iron line we concluded the spatial distribution of hot plasma is strongly non- 
uniform that should be taken into account in analysis of protons propagation in the 
GC. From the Suzaku data we got independent estimates for the diffusion coefficient 
of subrelativistic protons in the GC, which was in the range 10^^ — 10^^ cm^s~^ 
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1. Introduction 

This is a final paper of the series (Cheng et al. 2006; Cheng et al. 2007; Dogiel et al. 
2008; Dogiel et al. 2009abc) on energetic processes in the Galactic center (GC). In these papers 
we presented our interpretation of X-ray and gamma-ray emission from the Galactic center. 
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We supposed that all these phenomena had common origin and were initiated by accretion 
processes onto the central supermassive black hole. 

We showed that the the observed X-ray continuum and line emission from the GC might 
be produced by a flux of subrelativistic protons which resulted from an unbounded part of stars 
accreted onto the central black hole. We estimated the average energy of escaping protons to 
be about 100 MeV in order to produce a flux of hard X-ray emission in the energy range 
14-40 keV as observed by Suzaku/HXD (Yuasa et al. 2008). 

The quasi-stationary production rate of subrelativistic protons hardly exceeds Q ~ 2 x 
2q45 p];otons s~^ for the frequency of star capture ~ 10~^ — 10~^ years"^ (Syer & Ulmer 1999; 
Donley et al. 2002) and a fraction of escaped matter equaled ~ 50% of star masses (Ayal et al. 
2000). 

The unbounded fraction of stars escapes with subrelativistic velocities which correspond 
energies about 100 MeV for protons and 50 keV for electrons (see Dogiel et al. 2009a). The 
numbers of protons and electrons equal to each other. Therefore, fluxes of bremsstrahlung hard 
X-ray emission produced by protons and by electrons equal to each other too. But the lifetime 
of electrons with energies 50 keV is in about five orders of magnitude smaller than that of 
100 MeV protons and the electron bremsstrahlung flux is significant for a short time just after 
a capture event. Since it is supposed that relatively long time has passed after the last star 
capture, electron bremsstrahlung radiation is negligible. 

Our goal is to estimate the spatial diffusion coefficient, of subrelativistic particles 
near the GC, whose value is unknown. Below we derive its value from the spatial distribution 
of hard X-ray emission near the GC as observed by Suzaku. 

In section 2, we summarize an energetics of the Galactic Center diffuse emission. In 
section 3, we carry forward our model using non-uniform target gas distribution. Section 4 is 
devoted to explain how the model and the HXD data can be compared. In section 5 and 6, 
spatial distribution of 6.7 keV Fe line emission and hard X-ray continuum are compared with 
ones expected from the present model, deriving a confinement on a diffusion coefficient for the 
sub-relativistic protons. 

2. Components of Hard X-Ray Emission from the GC in the range 14 — 40 keV 

Koyama et al. (2007b) reported that the hard X-ray spectrum in the range 2 — 10 keV 
from the GC is naturally explained by a 6.5 keV-temperature plasma plus a power-law compo- 
nent with the photon index of F = 1.4. Latter on, Yuasa et al. (2008) found a prominent hard 
X-ray emission in the range from 14 to 40 keV whose spectrum is a power law with the spectral 
index ranging from 1.8 to 2.5. 

Different processes may contribute to the total emission from the GC. 

Thermal emission. The emissivity of thermal bremsstrahlung of a hot plasma e can be 
estimated from the equation (see e.g. Ginzburg 1989) 
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where the temperature is in Kelvin degrees. 

As follows from the analysis of Suzaku data the plasma temperature in the GC is 
constant and equals 6.5 keV (Koyama et al. 2007b). Then the thermal emissivity in the GC is 



Dogiel et al. (2009b) presented results of calculation for the 0.55° x 0.55° central region 
and showed that for the plasma density n = 0.2 cm~'^ and the temperature 6.5 keV the total 
thermal flux from this region in the energy rage 14 — 40 keV was about Fth — 2 x 10^^ erg s~^. 
The total non-thermal flux was estimated by the value Fnth — 3 x 10^^ erg s~^. 

Non-thermal component. Koyama et al. (2009) found that a combination of the Fe XXV- 
Ka and Fe I-Ka fluxes was proportional to the non-thermal continuum flux in the 5-10 keV 
energy band. They concluded that the total non-thermal flux consists of the two components, 
one of which is proportional to the intensity of 6.7 keV line, and, therefore, is produced in 
regions of hot plasma, while the intensity of the second component is proportional to that of 
the 6.4 keV line and is generated in regions of molecular gas. Warwick et al. (2006) reported a 
similar conclusion obtained from the XMM-Newton observations. 

In the model of X-ray production by subrelativistic protons this relation is naturally 
explained (see Dogiel et al. 2009bc). Continuum and line emission is produced in this case by 
interactions of subrelativistic protons with the background plasma and the molecular hydrogen 
in the GC. 

In spite of relatively small radius (r ~ 200 pc) the inner Galactic region contains about 
10% of the Galaxy's molecular mass, ~ ~ (7 — 9) x 10^ Mq. Most of the molecular gas is 
contained in very compact clouds (see the review of Mezger et al. 1996). Since the total mass 
of the molecular gas is much larger than that of the hot plasma, Mpi ~ 4 x 10^ Mq, one may 
assume that most of inverse bremsstrahlung flux produced by protons would come from regions 
filled by the molecular gas. However, as Dogiel et al. (2009c) showed, the mean free path of 
protons in dense molecular clouds is very short since the diffusion coefficient inside molecular 
clouds is quite small. Therefore, only a small fraction of the molecular hydrogen is involved 
into processes of inverse bremsstrahlung radiation. Thus, for the molecular cloud Sgr B2 with 
the mass about 10^ Mq it follows from calculations of Dogiel et al. (2009c) and the Suzaku 
data (see Koyama et al. 2007a) that the total bremsstrahlung flux in the range 2 — 10 keV is 
about -F2-10 — 10'^^ erg s^^. The total flux of the 6.4 keV line emission from Sgr B2 is about 
^6.4 ~ 1.2 X 10'^^ erg s^^ that gives the ratio -F2-io/-^6.4 ~ 8. For higher energies of X-ray photons 
our calculations show that the ratio of -Fi4_4o/-Fe 4 ~ 6.4. 

The total flux of the 6.4 keV line from the GC can be estimated from the ASCA data 
(Maeda 1998) and for the region 0.55° x 0.55° it is about ~ 1.3 x 10'^^ erg s~^, then we expect 
that the total non-thermal bremsstrahlung flux from proton interaction with the molecular gas 



Ane ~ 1.27 X 10 ^^n^(r) erg cm~^s~ 
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is about 10^^ erg s ^, and 70% of the total non-thermal bremsstrahlung emission in the range 
14 — 40 keV, ~ 2 X 10'^^ erg s"^, is generated in regions filled with the hot plasma. 

3. Qualitative effect of the plasma density variations 

In this section we present a qualitative effect of plasma density variations on the proton 
distribution in the GC. We remind that the spatial distribution of subrelativistic protons is 
described by the quasi-stationary diffusion equation 
d 
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where is the density of cosmic rays, D is the spatial diffusion coefficient, dE/dt = b{E) is the 
rate of energy losses which for subrelativistic protons is determined by the Coulomb collisions 
and has the form 

fdE\ _ 47rne^lnAi ^ k 

\dt).~ ~71' 
where Vp is the proton velocity, InAi is the Coulomb logarithm, and n is the plasma density, 
which is a function of the radius r in general case. The source function is taken in the form 

Q{E)=Q6{E-E,,,)6{r). (5) 

where Q ~ 3 x 10"^^ protons s~^ and Ecsc — 100 MeV is the injection energy of protons. 
For n = no = const, the spherically symmetric solution is 
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In the case when the plasma density is spatially variable as 
n{r) = no (^^^ , 

where tq is a scale parameter, the analytical solution for N{r,E) is 
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In figure 1 we show the difference between uniform for Uq = 0.2 cm~^ (solid line) and non- 
uniform cases for the parameters no = 10 cm~'^ (dashed-dotted line) and no = 3 cm~^ (dashed 
line). The scale parameter is taken to be ro = 25 pc. One can see that variations of the gas 
density change significantly the proton distribution. 
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Fig. 1. The spatial distribution of 70 McV protons in the GC for the case of uniform distribution (solid 
line) and that of non-uniform distribution equation (9) for uq ~ 10 cm^'^ (dashed-dotted line) and no = 3 
cm~'^ (dashed line). 

4. How to compare the model calculation and the observed HXD data? 

In section 5 and 6, we will compare spectra expected from the present model with the 
HXD spectrum observed around the GC (Yuasa et al. 2008). Before describing the comparison 
result, we present methods and assumptions which were used in the calculation. 

The intensity / of inverse bremsstrahlung emission of protons in any direction s is cal- 
culated from 




where the integration is along the line of sight s and Ey^ is the energy of photons. Here N{E,r) 
is the density of subrelativistic protons with the energy E at the radial distance r from the 
GC, V is the velocity of protons with the energy E, dai\^/dE^ is the cross-section of inverse 
bremsstrahlung radiation, and n(r) is the plasma density distribution in the GC. 

To compare such results of our model calculations, which are given in photons 
c'm~^keV~^s^^sr~^ with the Suzaku data of Yuasa et al. (2008) given in counts we should 
convolve the HXD energy and angular responses (see Mitsuda et al. 2007; Takahashi et al. 
2007; Kokubun et al. 2007) to the model intensity. The procedure is as follows: 

r /-40kcV 

/i4-4o=/ dn dEA{i,b)S{E)I{EJ,b), (12) 

Je,b J14keV 

where / is in {counts s~^) denotes the expected count rate in the direction determined by 
the galactic coordinates {£,b), Q is in (sr), I{E,£,b) in {ph cm~'^keV~^s~^sr~^) is the contin- 
uum intensity in the direction {i,b), S{E) and A{i, b) represent the effective area and angular 
transmission of the HXD, respectively. 
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The HXD data can still include a contaminating flux from a number of X-ray point 
sources besides the diffuse emission as noted in Dogiel et al. (2009b). The flux level of the 
contamination is calculated to be on the order of 10% by integrating the known LogN — LogS 
curve for X-ray point sources around the GC obtained with Chandra (Muno et al. 2009) over 
the luminosity of of 2 x 10'^^ — 1 x 10^^ erg s~^ in which range has been actually measured so 
far (in the "field" region of Muno et al. 2009). Although the contribution can be larger if the 
LogN — LogS curve is measured further smaller luminosities, in the present comparison with the 
model calculation, we simply ignore their contribution. The effect of this neglection is discussed 
later in section 6. 

5. Spatial Variations of the Gas Density in the GC 

The spatial distribution of the hot plasma in the GC can be derived from the distribution 
of the 6.7 keV iron line which traces the hot plasma. 

The origin of the K-a iron line from the hot plasma of the GC is pure thermal (see 
Dogiel et al. 2009b). Then for the surface brightness distribution of the 6.7 keV emission line 
in any direction s to the GC we have 

/6.7(s) oc J n'^{r)ds (13) 

s 

Maeda (1998) obtained a surface brightness distribution of the 6.7 keV emission line, 
/g.r, from the ASCA data as a function of the angle from the Galactic center, 6, 

/6.7(^,6) = /iexp(-M) + /2exp(-M) (14) 

UJl UJ2 

COS 9 = cos £ cos 6 

where 6 is the angle from {l,b) = (0,0), /i = 19.7 and J2=1.6 in photons crrr'^ sr~^ units, and 
oji = 0.42 and ci;2 = 15 in degrees. 

This distribution is shown in figure 2 by the thick solid line. 

From the 6.7 line distribution of Maeda (1998) we try to derive the plasma density 
distribution in the GC. Then we will use this density distribution for calculations of hard X-ray 
emission from the GC. We represent the hot plasma distribution by analytical functions as 

1 / r 

2 l^. 

and our goal is to find suitable parameters tq and a which give good correspondence of the 
observed ASCA data and that obtained with the density distribution (15). 

In order to derive the model spatial variations of the 6.7 keV line emission along the 
Galactic plane we should integrate along the direction s shown in figure 3. The angle between 
the direction to the GC and the line of view is the Galactic longitude i. The circle line around 
the Galactic center defines a sphere filled with a hot plasma whose radius is supposed to be 
about a = 200 pc in accordance with the X-ray observations (see, e.g. Koyama et al. 1996). 
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Fig. 2. The latitude distribution of 6.7 keV line as from the ASCA data byMaeda (1998) is shown by the 
solid line. Model spatial distribution of the surface brightness profile of the 6.7 keV emission line expected 
from equation (17) for plasma distributions: ro = 25 pc and a = 0.7 (dashed-dotted line), and uniform 
plasma distribution (thin solid line). 
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Fig. 3. The schematic view of the GC from Earth. The arrow line shows the distance between Earth and 
the GC, the solid line is the line of view, the solid circle shows the region of hot plasma around the GC; 
the two radius- vectors mark positions of ipmax and ipmin = tt — ■(/'max- 

The radius-vector r from the GC along the hne of sight as a function of the Galactic longitude 
i and the angle between the line coinciding Earth an the GC (arrow line) is described as 

sin£ 



r{i,tlj) = d- 



cos ip 



(16) 



Then with the density distribution (15) the integration along the line of sight in the direction 
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and the angles i^/'max shown by the two radius-vectors in figure 3 is 



^max = arccos 
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(19) 



Here d = 8 kpc is the distance between Earth and the GC which equals 8 kpc (Reid 1993). From 
equation (17) and the ASCA 6.7 keV spatial variations obtained by Maeda (1998) we derived 
parameters vq and a, that gives ro = 25 — 75 pc and a = 0.7— 1.7. The ASCA spatial variations 
of 6.7 keV line are shown in figure 2 by the thick solid line. As an example we show also in 
this figure the model spatial variations of the 6.7 keV line intensity derived for tq = 25 pc and 
a = 0.7 (dashed-dotted line). The coincidence between calculations and data are reasonably 
good. For comparison we showed also in this figure the expected 6.7 keV line spatial variation 
for the case of uniform plasma distribution (thin solid line). One can see that the ASCA data 
can hardly be described by a model with the uniform plasma distribution. 

The parameter rio of equation (15) can be roughly estimated by comparing the model 
emission (6.5 keV thermal + non-thermal) with the observed spectrum of the hard X-ray 
emission around the GC (Yuasa et al. 2008). To perform such a comparison, we convolved 
the model spectrum with the HXD response as described in section 4, and then, overlaid it 
over the HXD spectrum by adjusting the model normalization factor, uq in equation (15), so 
that the model reproduces the HXD fiuxes. Figure 4(a) shows a result of the convolution, and 
the model well reproduce the observed spectrum in the 14 — 40 keV band. Although we only 
present the result for one parameter set explained in the caption of figure 4, the model spectral 
shape does not strongly depend on the parameters, and within possible parameter ranges, no 
was calculated to be ~ 0.1 — 0.3 cm~'^. 

As one can see from these figures the contributions of thermal and non-thermal emission 
into the total fiux from the GC are almost equal to each other in this energy range (their ratio 
is 2:3). 

6. Diffusion CoefHcient of Cosmic Rays in the Galactic Center 

In order to calculate the spatial distribution of subrelativistic protons in the GC we 
should take into account energy losses in regions of hot plasma and in the regions of neutral 
gas surrounding the GC whose average density was taken to be Uh = 1 cm""^. 

As follows from the analysis of Yuasa et al. (2008) the spatial distribution of the intensi- 
ties of thermal and non-thermal emission in the GC are similar to each other. Since the first one 
is proportional Ith oc / n^{r)ds and the second one is proportional to /nth oc / n{r)Np{r)ds, where 

r r 

Np{r) is the spatial distribution of subrelativistic protons, it means that the spatial distribution 
of subrelativistic protons does not differ strongly from that of plasma. 

The expected spatial distributions of the 14 — 40 keV X-ray emission in the GC derived 
from equation (11) for = 3 x lO^^cm^s"^ (solid line) and for D = 3 x lO^^cm^s"^ (dashed line) 
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Fig. 4. (a)Thc 14 — 40 kcV X-ray spectrum observed by HXD/PIN around the GC (taken from Yuasa 
ct al. 2008). The model spectra for 6.5 keV thermal emission, inverse bremsstrahlung, and sum of them 
are also shown in dashed, dashed-dotted, and solid line. The model spectra were calculated using such 
parameters as ro = 50 pc, a = 1.1, D = lx 10^^ cm^/s, and no = 0.13 cm~'^. (b)The same model spectra as 
(a) are shown without convolving the HXD /PIN energy response but with integration over the PIN FOV 
(0.°55 X 0.°55, effectively). 

and the Suzaku data obtained by Yuasa et al. (2008) are shown in figure 5 ^. Our calculations 
show that spatial distribution of the non-thermal 14 — 40 keV X-ray emission does not diff'er 
significantly from each other. However, the necessary proton production rate Q is an essential 
function of the diffusion coefficient. Thus, for D = 3 x lO^^cm^s"^ the value of Q needed to 
reproduce the Suzaku data equals lO^^protons s~^ while for D = 3x lO^^cm^s"^ the production 
rate is necessary to be Q = 2.3 x lO^^protons s"^. We remind that the average rate of proton 
production by accretion on the black hole is no more than Q = 3 x lO^^protons s~^. Since the 
value of Q is restricted, then permitted values of D are confined within a relatively narrow 
range around D ^ 10^^ cm^ s~^. 

However, these estimates were obtained in assumption that almost 100% of the 14 — 
40 keV X-ray fiux from the GC is generated by subrelativistic protons. As follows from, e.g., 
Warwick et al. (2006); Revnivtsev et al. (2009) a significant part of the GC hard X-ray fiux 
may be due to faint point sources. Then a smaller production rate of protons in the GC is 
required, and, hence, higher values of the diffusion coefficient cannot be excluded. 

^ In figure 5, the data point at I = 0.°55 is largely discrepant from the model, and even from the other data 
points. As noted in Yuasa et al. (2008), this deviation could be due to underestimation of contaminating 
signals from known bright point source inside the HXD/PIN FOV. Therefore, this point should be ignored 
in the present comparison between observed data and the model calculation. 
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Fig. 5. The spatial distribution of inverse brenisstrahlung flux in the energy range 12 — 40 kcV averaged 
over the HXD field of view with the Suzaku data. Solid line: D = 3x lO^^cm^s"^, Q ~ lO^^protons s~^. 
Dashed line: D = 3x lO^^cm^s"^, Q = 2.3 x lO^^protons s~^. The Suzaku data were taken from Yuasa 
et al. (2008) 

7. Discussion and Conclusion 

Analysis of X-ray data (Koyama et al. 2009) demonstrated that continuum non-thermal 
emission in the GC consisted of two components one of which was proportional to the intensity 
of 6.7 keV iron line which traced distribution of hot 6.5 keV plasma, and the other one was 
proportional to the intensity of 6.4 keV plasma which traced distribution of cold molecular 
hydrogen. If this non-thermal emission is due to inverse bremsstrahlung emission of subrela- 
tivistic protons, then their interactions with hot and cold fractions of the interstellar medium 
are equally important. Our estimation show that about 30% of the total non-thermal flux from 
the GC in the range 14 — 40 keV is generated in regions of cold gas while the rest is produced 
by proton interaction with hot plasma. The reason is that only a small fraction of molecular 
gas is involved into processes of bremsstrahlung emission because of short mean free path of 
protons in molecular clouds. 

From the X-ray data obtained with ASCA (Maeda 1998) it was concluded that the 
plasma density distribution in the GC is strongly non-uniform, e.g. the model distribution of 
6.7 keV line in the GC obtained with the plasma profile 



is very close to that from ASCA observation. This effect of density variations should be included 
into analysis of protons propagation in the GC. 

From the total spectrum of X-ray emission in the range 14 — 40 keV from the 0.55°- 




(20) 
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diameter central region and the spatial variation of this emission within 1.5°-diameter central 
region as observed by Suzaku and ASCA we derived the characteristic diffusion coefficient of 
subrelativistic protons. It was shown that the diffusion coefficient of subrelativistic particles 
is about the value of 3- 10^^ cm^s~^. For higher values of the diffusion coefficient a higher 
production rate of subrelativistic protons is necessary which cannot be provided by accretion. 
However, this estimation of D one should consider as a lower limit. The point is that we assumed 
here that all hard X-ray emission in the 14 — 40 keV band is produced by subrelativistic protons 
in the hot fraction of the interstellar gas. However as stated above a part of this emission may 
be produced in molecular clouds. Besides, we do not know exactly which part of this emission 
is due to unresolved point sources (see discussion in Dogiel et al. 2009b). Therefore, we cannot 
exclude that the diffusion coefficient in the GC is close to its average value in the Galaxy, i.e. 
D ~ 10^'^ cm^ s~^. Thus, we determine the range of values of the diffusion coefficient in the GC 
as 10^6 - lO^^cm^s-^ 
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